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We present the first interpretation of the new isotropic diffuse 7-ray background (IGRB), measured by the 
Fermi Large Area Telescope (LAT), based on a statistical analysis. We demonstrate that the 7-ray emission 
from unresolved active galactic nuclei and star forming galaxies is consistent with the Ferrm-LAT IGRB data 
within the uncertainties both on the choice of the Galactic diffuse emission model and on the 7-ray emission 
mechanism of these sources. Furthermore, adding to the extragalactic sources the contribution from a smooth 
Galactic halo of annihilating weakly interacting dark matter (DM) particles, we are able to set stringent limits 
on the DM annihilation cross section. Finally, we demonstrate that the addition of DM can significantly improve 
the fit to IGRB data. 


1. Introduction 

Recently, the Large Area Telescope (LAT) on board 
the Fermi Gamma-ray Space Telescope (Fermi) has 
published a new measurement of the isotropic diffuse 
7-ray background (IGRB) and the extragalactic 7-ray 
background (EGB) in the energy range 100 MeV-820 
GeV Q]. The origin of this 7-ray residual represents 
one of the most mysterious open problems in astro¬ 
physics. The IGRB is usually associated to the 7-ray 
emission from unresolved, namely not detected by the 
Fermi- LAT, blazars, misaligned active galactic nuclei 
(MAGN) and star forming (SF) galaxies [21 IS HI- The 
most powerful Galactic contributors of the IGRB are 
expected to be pulsars, due to the large sample of 
detected sources. However, very recently in Ref. [5j 
the 7-ray emission from high-latitude (|6| > 20°) un¬ 
resolved pulsars has been derived to account for, at 
most, the 1% of the IGRB. 

The annihilation of dark matter (DM) particles, in the 
halo of the Milky Way (MW) and in external galax¬ 
ies, constitutes a possible exotic mechanism for the 
production of 7 rays. Indeed, one of the most promis¬ 
ing strategies for the search of DM, in the scenario of 
weakly interacting massive particles (WIMPs), is the 
indirect detection through 7 rays produced from the 
annihilation of WIMPs. 

We first explore in Sec. [2] at which extent the astro- 
physical populations may explain the EGB and IGRB 
data and then in Sec. [3] we add the contribution from 
a DM Galactic halo to constrain the DM annihila¬ 
tion cross section. We use a statistical fitting proce¬ 
dure which includes both the Fermi-LAT data errors 
and the theoretical uncertainties on the 7-ray emis¬ 
sion from astrophysical sources. For the first time we 
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show the effect of the choice of the Galactic diffuse 
emission (GDE) model, used to derive the IGRB and 
EGB data, and how this affects the results on DM. 

2. Fermi-LPJ IGRB data explained with 
astrophysical sources 

At least the 10% of the 7-ray photons, detected 
at latitude \b\ > 20° by the Fermi- LAT, are emit¬ 
ted from Galactic and extragalactic resolved sources. 
Indeed, the Fermi- LAT catalogs contain thousands of 
point sources. One of the most natural explanation 
for the origin of the IGRB, is that it arises from the 
superposition of a numerous population of unresolved 
sources with a flux lower than the point source sensi¬ 
tivity threshold of the LAT [21 [5] . 

A diffuse 7-ray flux has been predicted for various 
source populations. A large fraction of the IGRB is 
expected to come from active galactic nuclei (AGN) 
which are divided, according to the orientation of the 
jets with respect to the line of sight, into blazars and 
MAGN. Blazars are the most numerous population 
in the Fermi- LAT catalogs and are divided into BL 
Lacertae (BL Lac) objects and flat-spectrum radio 
quasars (FSRQs) according to the absence or pres¬ 
ence of strong broad emission lines in their optical/UV 
spectrum, respectively. The 7-ray emission from un¬ 
resolved blazars has been estimated to be around 20- 
30% of the integrated IGRB above 100 MeV (see e.g. 
mm®- MAGN are only a dozen in the 7-ray cata¬ 
logs but the unresolved counterparts are expected to 
be much more numerous than blazars. The diffuse 
7-ray emission from MAGN is expected to be aroud 
20-30% of the integrated IGRB above 100 MeV OS]. 
Finally, SF galaxies are predicted to have a numerous 
unresolved population and their contribution has been 
derived to be from a few % up to almost the totality of 
the IGRB (see e.g. [31 HO])- In particular the Fermi- 
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LAT Collaboration [3] has deduced that the contribu¬ 
tion of unresolved SF galaxies is among 4% and 24% 
of the integrated IGRB above 100 MeV and they have 
considered two different Spectral Energy Distribution 
(SED): a power law (PL) and a Milky Way (MW) 
function. 

Truly diffuse processes, like ultra-high energy cosmic- 
ray interactions with the extragalactic background 
light or intergalactic shocks, are other possible 7 - 
ray emission mechanisms (see Ref. m and references 
therein). However, there are still quite large uncer¬ 
tainties associated to these processes and a subdom¬ 
inant contribution to the IGRB is theoretically pos¬ 
sible (see e.g. Refs. P2] na IT4]h Therefore, we do 
not take into account these contributions in the rest 
of the paper. 

We consider in our analysis the 7 -ray emission from 
unresolved SF galaxy, MAGN, BL Lac, and FSRQ 
populations as derived in In Fig. [T] we 

display the Fermi -LAT IGRB and EGB data derived 
with the Model A of GDE [T] , together with the 7 -ray 
fluxes predicted for AGN and SF galaxies. It is evident 
that both the IGRB and EGB data are consistent with 
the emission from these extragalactic populations. 

We now determine with a chi-square statistical 
method at which extent AGN and SF galaxies can 
explain the IGRB and EGB data. We perform a fit 
to the IGRB and EGB measurements with a \ 2 func¬ 
tion including both the errors of the Ferrm-LAT data 
and the theoretical uncertainties of the 7 -ray emis¬ 
sion from astrophysical sources (see for all the details 
ini)- The theoretical uncertainties on AGN mainly 
produce a renormalization of their average unresolved 
7 -ray flux while the SF galaxy contribution contains 
also a large uncertainty due to the SED. We thus con¬ 
sider, at the first order, the l-cr error on the AGN 
and SF galaxy unresolved emission as the width of 
the bands shown in Fig. [l] while the average fluxes 
are represented by the curves in Fig. |T] Both the SF 
galaxy SED, namely the MW and the PL models, are 
taken into account to include also our ignorance on 
the spectral shape of these sources. 

We explore the possibility that the theoretical predic¬ 
tions adopted may be affected by an additional un¬ 
certainty on the spectral shape. This possibility is in¬ 
cluded by varying the power-law index of AGN taking 
into account the relevant l-cr error on this parameter 
(see for further details [Hj). The uncertainties on the 
SF galaxies SED is already considered by performing 
all our analysis with both the MW and PL models. 
The results are shown in Tab. [Tl where we display 
the reduced chi-square (% 2 = x^/d.o.f.) for the fit 
on the IGRB and EGB data. The analysis is per¬ 
formed on the Fermi -LAT data derived with each of 
the benchmark GDE models considered in [T] and la¬ 
belled with A, B and C. The 7 -ray emissions from 
unresolved BL Lacs, FSRQs, MAGN and SF galaxies 
are able to explain the high-latitude IGRB and EGB 


X 2 

IGRB (MW) 

EGB (MW) 

IGRB (PL) 

EGB (PL) 

MODEL A 

MODEL B 

MODEL C 

1.72; 1.56 

1.33; 1.32 

0.82; 0.84 

0.95; 1.02 
1.57; 1.72 
0.60; 0.60 

3.20; 2.54 
2.30; 1.96 
1.67; 0.95 

1.41; 1.36 
1.83; 2.06 
0.77; 0.84 


Table I The reduced chi-square value y 2 = x 2 /d.o.f. for 
the fits performed using only the normalizations as free 
parameters (left numbers in each column) or varying also 
the slope of the spectra (right numbers in each column). 
The analysis is performed with both the MW and PL SF 
galaxy models and for the three Galactic foreground 
models of the IGRB and EGB data p]. 


data with no need for significant adjustments of the 
average parameters. A better agreement is provided 
by the MW modeling of the SF galaxy emission and 
with the Model C of the GDE [T] . Indeed, the PL 
model of SF galaxies is in some tension with the data 
sets obtained with GDE Model A and B (see Tab. [I]). 
The choice of the Galactic foreground model has a 
large relevance on the goodness of the fit. We display 
in Fig. [2] the fluxes corresponding to the best fit con¬ 
figuration with the IGRB and the EGB data obtained 
with Model A of the GDE. It is remarkable that each 
contribution has a very different shape but they add 
in a way that their sum is consistent with a power-law 
with an exponential cutoff. 


3. Constraints on a DM contribution to 
the Fermi- LAT IGRB data 

A possible contribution to the high latitude 7 - 
ray IGRB could arise from annihilating DM particles 
present both in the halo of the MW and in external ha¬ 
los m- DM can produce 7 rays both directly (the so- 
called prompt emission), or indirectly via the inverse 
Compton scattering of the electrons and positrons, 
produced by the DM annihilation, off the ambient 
light of the interstellar radiation field. In order to 
simplify the discussion we consider only the DM dis¬ 
tributed in a halo of the MW. Moreover, we do not 
consider any specific particle physics model and we fix 
the branching ratio equal to 1 for any of the discussed 
annihilation channels. The photon and electron spec¬ 
tra have been calculated using the Pythia Montecarlo 
code for DM annihilations into e + e _ , n + , t + t~, bb , 
tt and W + W~ channels. We have assumed an Einasto 
DM profile with a local DM density of p© = 0.4 
GeV/cm 3 . For all the details about the 7 -ray flux 
due to DM annihilation in the halo of the MW we 
refer to mm- 

As a first analysis, we derive upper limits for the 
DM annihilation cross section ( av ), fitting the IGRB 
and EGB data with the astrophysical sources dis¬ 
cussed in Sect. [2] and the addition of a Galactic DM 
halo. The results are shown in Fig.[3]for the DM anni- 
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IGRB composition with MW SF model 


EGB composition with MW SF model 
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Figure 1: Left (right) panels show the 7 -ray emission from unresolved (unresolved-l-resolved) sources, together with the 
IGRB (EGB) data HQ. Lines and relevant uncertainty bands represent the contribution from the following source 
populations: MAGN (orange dashed), BL Lacs (green dotted), FSRQs (grey double dot-dashed), MW model of SF 
galaxies (purple dot-dashed), and the sum of all the contributions (blue solid). IGRB and EGB data have been derived 
with Galactic foreground Model A. 


FIT TO FERMI-LAT IGRB MODEL A 


FIT TO FERMI-LAT EGB MODEL A 




Figure 2: The 7 -ray flux best fit on Model A of IGRB (left panels) and EGB (right panels) data is shown for: BL Lac 
(dotted green), MAGN (dashed brown), SF (dot-dashed purple), FSRQ population (dot-dot-dashed black), the sum of 
AGN and SF (solid blue). We display with a cyan band also the uncertainty associated to the Fermi -LAT GDE model 

0 - 


hilation channels listed above, different confidence lev¬ 
els (C.L.s) and considering the Model A of the Galac¬ 
tic foreground. In the case of the bb DM annihilation 
channel and for masses lighter than 30 GeV, the upper 
bounds on ( av) are below the thermal relic value while 
at 10 TeV our analysis excludes (av) > lCV 24 cm 3 /s. 
Moreover, the limits obtained for the DM annihilation 
into t + t~ are quite stringent: the thermal relic cross 
section is excluded at 3-cr C.L. up to a DM mass of 
about 330 GeV, while at m x ~ 1 TeV the bound is 
around 10 _25 cm 3 /s. As expected, the upper bounds 
obtained with a fit to the EGB data are very similar 
to the ones obtained with the IGRB data and in gen¬ 
eral are only slightly looser. 

For the first time we explore the impact of the GDE 
model in the estimation of upper limits on the DM 
annihilation cross section. The results are illustrated, 
for a 2 -a C.L and for the annihilation channels e + e~ 
and bb, in Fig. [ 4 ] The choice of the GDE model turns 
out to have a significant role in the values of the up¬ 
per limits on (av ). The results vary on average within 
a factor of two but the differences can reach a factor 
of 10 in the case of bb at about m x ~ 200 GeV. The 
bounds obtained on the IGRB with Model B of the 


GDE are always looser with respect to the ones de¬ 
rived with Model A or C. 

As a second analysis, we attempt to identify DM 
configurations which can significantly improve the fit 
to the IGRB data. In this part we perform a fit to 
the Fermi -LAT data with extragalactic sources and a 
DM component as done before, but letting the WIMP 
DM mass m x and (av) varying simultaneously. Our 
results are displayed in Fig. [HJ for the bb DM annihi¬ 
lation channel and for the IGRB data associated to 
the three Galactic foreground models. In the case of 
Model A and B, we obtain closed regions up to 3- 
a C.L. with the best fits located around m x ~ 5-20 
GeV and (av) ~ 1-3 ■ 10 _26 cm 3 /s while for Model C 
the 1 -a C.L. closed region opens up at already 2-a 
C.L., translating the results into upper limits. The 
addition of a DM component is almost irrelevant for 
IGRB Model C while, in the case of leptonic or bb 
channels, improves the IGRB fit, for models A and B, 
with Ay 2 = Xastro Xastro+DM > 6-5, where Xastro and 
Xastro+DM are the best fit x 2 with only astrophysical 
sources or with the addition of a DM component. In 
the case of e + e _ , t + t~ and bb channels, the best fit 
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Figure 3: Upper limits on (crv) as a function of the DM mass m x for bb (top panels) and t + t~ (botton panels) 
channels. The upper limits are derived with a fit to the IGRB (left panels) and EGB (right panels) data, within GDE 
Model A. The 3-cr, 2-er and 1-cr C.L.s are shown with solid, dashed and dotted curves, respectively. The horizontal 
dotted line specifies the thermal relic annihilation cross section value. 


e + e" 2-o ULs changing Galactic foreground 


bb 2-o ULs changing Galactic foreground 




Figure 4: 2 -a C.L. upper limits on the annihilation cross section for the e + e (left panels) and bb (right panels) 
channels obtained with the IGRB data derived with three different GDE models [T]. 


values of the DM mass range from few GeV up to 20 
GeV and the annihilation cross section values are close 
to the thermal one while the channel requires 

(crv) ~ 1-3 • 10 _25 cm 3 /s. It is remarkable that includ¬ 
ing a DM component does not require the standard 
astrophysical contributions to differ significantly from 
the average emission m- Therefore, a DM compo¬ 
nent can very well fit the IGRB data with a realistic 
unresolved emission from extragalactic sources. 

The results illustrated in Fig. [5] demonstrate that a 
DM contribution to IGRB may significantly improve 


the fit to the IGRB with respect to the interpreta¬ 
tion with only astrophysical source populations. How¬ 
ever, the significance of this potential exotic signal 
strongly depends on the choice of the Galactic fore¬ 
ground model considered to derive the IGRB spec¬ 
trum. It is evident in Fig. [5] that for Model A and B 
of the GDE we may have an hint of DM up to 3-cr C.L. 
while for Model C for 2 and 3-<r C.L. we can only set 
upper limits on the annihilation cross section. These 
results confirm how a deep knowledge of the GDE is 
mandatory to unveil a DM contribution in the IGRB 
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Figure 5: Contour plots in the (av)-m x plane for bb channel calculated with a fit of the astrophysical sources and a DM 
component on the IGRB data. The black dots refer to the values of the best fit, and the (closed or open) green, blue 
and pink regions indicate the 1, 2 and 3-<r C.L. Left, center and right panels correspond to the case of Model A, B and 


C of the IGRB. 


data. In Fig. [ 6 ] we illustrate the best fit configuration 
for the bb DM annihilation channel and Moled A of the 
GDE. This is the case with the largest significance of 
DM with respect to the best fit with only astrophysi¬ 
cal populations with Ay 2 = 18.9. In this specific best 
fit, the DM mass is m x = 8.2 GeV, and the annihila¬ 
tion cross section is (av) = 1.4 • 10” 26 cm 3 /s. It is 
clear in Fig. [ 6 ] that this configuration of extragalactic 
sources and DM reproduce very well the IGRB and 
EGB data. 

The upper bounds reported in Figs. [3] and [4] improve 
the results derived in El. in the so-called ‘best-fit’ 
scenario, by a factor of ~3 at m x ~ 10 GeV and a 
factor of 30 at m x ~ 10 TeV. Our limits also improve 
significantly the upper limits on (av) derived by the 
Fermi -LAT analysis for a Galactic halo of DM [IS] 
and with the analysis of 25 dwarf Spheroidal galaxies 
[19] . We notice that for m x ~ 10 TeV the upper 
bounds found in [20] by the H.E.S.S. Collaboration, 
which are optimized at energies larger than about 1 
TeV, are of the same entity as ours for the leptonic 
channels, while for hadronic channels they are about 
one order of magnitude weaker. In addition, we derive 
similar results with respect to the very recent analysis 
performed by the Fermi -LAT Collaboration in Refs. 
PH [22]. 


4. Conclusions 

We have performed the first detailed statistical 
analysis for the interpretation of the recent IGRB 
data, measured by the Fermi -LAT Collaboration [Tj. 
We first test the hypothesis that a numerous sam¬ 
ple of unresolved extragalactic sources may explain 
the Fermi-LAT IGRB data. For the first time a \ 2 
function which includes the data errors and the the¬ 
oretical uncertainties on the 7 -ray emission from BL 
Lacs, FSRQs, MAGN and SF galaxies has been con¬ 
sidered. The theoretical uncertainty from each of this 
extragalactic population has been first parametrized 


only with a renormalization factor and then general¬ 
ized adding also a possible change in the photon in¬ 
dex of the 7 -ray SED. Since the spectral shape of SF 
galaxies is not well constrained we have considered 
both the MW and the PL models. In our results the 
IGRB and the EGB data are well fitted by the unre¬ 
solved emission from AGN and SF galaxies with best 
fit parameters close to the average theoretical values. 
We also demonstrate how the choice of the Galactic 
foreground model, used to derive the IGRB and EGB 
data, affects the results. 

We explore also a possible contribution from the an¬ 
nihilation of WIMP DM particles distributed in the 
halo of our Galaxy. As a first analysis, we derive up¬ 
per limits on the DM annihilation cross section, com¬ 
bining the 7 -ray emission from astrophysical sources 
and DM in order to fit the IGRB and EGB data. The 
upper bounds calculated with this method are strin¬ 
gent and rule out the thermal relic cross section for a 
wide range of DM mass values for the bb and t + t~ 
annihilation channels. 

We finally derive DM configurations which improve 
the fit to the IGRB data with respect to the case 
with only astrophysical sources. The best fit DM mass 
ranges from few GeV up to 20 GeV and the annihila¬ 
tion cross section (av) values are close to the thermal 
one. A DM component may fit, together with AGN 
and SF galaxies, very well the IGRB and EGB data 
with best fit parameters for the astrophysical popula¬ 
tions close to the average theoretical values. However 
depending on the Galactic foreground model and the 
value of the C.L., we obtain a possible hint of DM or 
we set only upper limits. 

Our results show how crucial is the IGRB in the study 
of the extragalactic source populations and for DM 
searches. It is today a powerful tool to constrain the 
DM properties and with a future improvement in the 
knowledge of the GDE and of the unresolved emission 
from AGN and SF galaxies, may probe a DM contri¬ 
bution to the 7 -ray sky. 
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IGRB best fit with DM, bb channel 


EGB best fit with DM, bb channel 




Figure 6 : The differential 7 -ray flux for the unresolved (left panels) and entire (right panels) BL Lac, FSRQ, MAGN, 
SF galaxy populations and the DM contribution with bb channel (splitted into the prompt and the ICS emission) as 
fixed by the best fit to the Model A of the IGRB and EGB data. 
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